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SUMMARY 

A new viscous- inviscid interaction procedure is presented 
which is applicable to separated flow’s. The new procedure is 
simple, converges rapidly, and does no* require numerical smooth- 
ing and underrelaxat ion , at least in the cases computed thus far. 
Calculations are presented for the low speed separated flow in 
the juncture region between an axisyimaetric body and sting. The 
viscous computation is done with an inverse boundary-layer pro- 
cedure which was previously developed. The inviscid computation 
is made with an axisysssetric transonic code called BAXRGD. The 
main advantage -of the new interaction procedure is that it com- 
bines an inverse boundary-layer technique, which is applicable 
to separated flows, with an existing inviscid analysis code with 
only a slight boundary condition change required in the inviscid 
code. 


INTRODUCTION 

An accurate prediction of the aerodynamic forces on a vehicle 
requires that the interact ior. ’.etaeeu the viscous and inviscid 
flow be taken into account, particularly near the aft end where 
the thickness of the viscous layer is significant and separation 
often occurs. If separation is present, conventional matching of 
the boundary-layer and inviscid flow is divergent since the 
boundary-layer solution is singular at the separation point, 
resulting in a displacement body with a vertical slope. This 
breakdown due to separation has been avoided in several conven- 
tional boundary- layer approaches by making various assumptions 
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about the separated region such as that presented in reference 1 
for supercritical airfoils and reference 2 for axisymmetric boat- 
tail flows. This singularity can be eliminated by solving the 
boundary-layer equations inversely; that is, the displacement 
thickness is prescribed and the pressure distribution is deduced 
from the boundary-layer solution. In reference 3, a recently- 
developed inverse boundary-layer formulation and finite-difference 
solution technique is presented and its applicability to turbulent 
separated flows is demonstrated. This new boundary-layer procedure 
can be used in either the direct mode (pressure prescribed, attached 
flow) or inverse mode (displacement thickness prescribed), and is no 
more complicated than conventional solution techniques. 

The present difficulty in using an inverse boundary- layer 
technique is the coupling of such a procedure to a simultaneous 
computation of the inviscid flow to iteratively solve for the 
viscous-inviscid interaction. Some progress has been made in this 
area with the development of the interaction procedures given in 
references 4 and 5 for low-speed laminar flows and references 6 
and 7 for transonic turbulent flows. However, significant improve- 
ments are needed in the first three of these techniques since they 
have various limitations and have not received widespread use. 

The technique developed by Thiede in reference 7 appears to be 
quite promising as it combines an inverse integral boundary- layer 
technique with a relaxation solution of the exact potential equa- 
tion through a novel matching procedure presented in reference 8. 
Unfortunately, Thiede only presented the outline of this inter- 
action procedure; further details are required in order to 
implement this scheme. 

The purpose of this report is to present a new interaction 
procedure which permits the inverse boundary-layer procedure re- 
ported recently in reference 3 to be combined with an existing 
inviscid flow analysis. This new interaction technique is demon- 
stra*ed by presenting solutions for the separated turbulent flow 
over axisymmetric body sting juncture which is schematically 



shown in figure 1. The inviscid calculations were made with an 
axisymmetric transonic code called EMBOD. This code, which 
numerically solves the exact potential equation by relaxation, 
was originally developed by South and Jameson (ref. 9) and later 
modified and documented in reference 10. These calculations are 
for low-speed flow only; however, work is underway to extend 
these calculations to transonic speeds using the compressible 
inverse boundary- layer procedure in reference 3. This new inter- 
action procedure, which is simple and converges rapidly, should 
be equally applicable to two-dimensional flows such as that over 
an airfoil, provided an appropriate inviscid code is used. 

SYMBOLS 

C f skin friction coefficient 

C p coefficient of pressure 

H boundary- layer shape factor 

i interaction iteration index 

Be^ Reynolds number based on free stream velocity and body length 
r Q body radius 

u velocity component parallel to body surface 

u e velocity component parallel to body surface at edge of 

boundary layer 

u e v u e Predicted by viscous calculation 

u e predicted by inviscid calculation 

I e 

v n normal velocity component at body surface 
x coordinate measured along body surface 

6 displacement thickness 

0 momentum thickness 

6) relaxation parameter 
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p density 

t . shear stress at body surface 

w 

VISCOUS ANALYSIS 

The inverse boundary- layer technique used in the present cal- 
culations is the 6 -prescribed procedure presented in refer- 
ence 3, which is a finite-difference technique for the two- 
dimensional incompressible boundary- layer equations. For the 
present application, the Mangier transformation (ref. 11) is 
used to convert the two-dimensional boundary-layer solution to 
the required axisymmetric solution. In the present investigation, 
the Crank-Nicolson finite-difference scheme used in reference 3 
was replaced with a first-order accurate scheme to suppress the 
streamwise oscillations which occurred in the body-sting juncture 
region. The first-order accurate scheme has greater streamwise 
damping than the second-order accurate Crank- Ticolsoa scheme, 
which has neutral stability in the wall region. This same modi- 
fication was used in a previous investigation reported in refer- 
ence 12 in which oscillations resulted when the Crank-Nicolson 
scheme was used to solve the boundary-layer equations subject 
to discontinuous wall cuotion. 

In the present calculations, even with the first-order scheme , 

some oscillations still occurred in the solution unless the term 

$ 

X 

-£■£-, which appears in the x -momentum equation, was differenced 

in its conservative (unexpanded) form. With m used to denote 
the x-grid point index, the conservative differencing of this 
term gives 


3 u6* 

dx 


m 


(u6*) m - (u6 *) ffl _i 
Ax 


+ 0( Ax) 


( 1 ) 
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where O(Ax) is the first-order truncation error. In the pre- 
vious study (ref. 3) this issue did not arise since both the 
conservative and nonconservative form of this term results in 
the same finite-difference expression, which is given as follows, 
when the Crank-Nicolson scheme is used: 


3u6* 

dx 


O-S 


(Ufi\ - 

Ax 


+ 0(Ax 2 ) 


( 2 ) 


INVISCID ANALYSIS 

The inviscid calculations were done with a program called 
RAXBOD (ref. 10) which solves the finite-difference form of the 
exact potential equation by successive line overrelaxation for 
the transonic (or subsonic) inviscid flow over axisymmetric 
bodies in free air. For bodies with a sting this program uses 
a body-normal coordinate on the forefoody up to the first hori- 
zontal tangent and a sheared cylindrical coordinate system aft 
of that point. This coordinate system seeas particularly well 
suited to bodies which have sharp corners as demonstrated by 
South and Jameson (ref. 9). A stretching is applied to both 
the normal and tangential coordinates such that the infinite 
physical space is mapped to a finite computational space. A 
typical mesh distribution is shown in figure 2. This program 
is easy to use and has been combined by several investigators 
such as Wilmoth (ref. 2) with direct (pressure prescribed) 
boundary-layer calculations. 

INTERACTION PROCEDURE 

The iterative procedure, which is shown in figure 3, begins 
with an initial guess for 6 , the displacement thickness distri- 
bution, which is then input to the inverse boundary-layer code; 
and the resulting solution yields the viscous edge velocity. The 
displacement effect of the viscous layer on the external inviscid 
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flow is represented by prescribing, in the inviscid flow computa- 
tion, an injection (or suction) velocity at the body surface which 
is given by 


A 

dx 


< 6 Vo> 


<3) 


for axisynsmetric incompressible flow. This representation of the 
boundary layer is equivalent, but preferable, to adding the dis- 
placement thickness to the original body; in the latter case a 
recomputation of the inviscid mesh and a complete restart of the 
inviscid solution is required for each interaction cycle. Use of 
equation (3) affects only the surface boundary condition in the 
inviscid calculation and thus does not require significant modi- 
fications be made to the inviscid code. The output from the 
inviscid calculation which is made with RAXBGD in the present 
case, is the inviscid surface tangential velocity. 

The next step in the interaction procedure is the computa- 
tion of the new displacement thickness distribution which, as 
shown in figure 3, is deduced from the mismatch of the viscous 
and inviscid tangential velocities by 


u e 

S new ” 5 old u 

e I 


(4) 


This update procedure for the displacement thickness is the key 
link which allows an inverse boundary-layer computation to be 
combined with a direct inviscid flow analysis. The iterative 
process continues until convergence is obtained when the new and 
old displacement thickness distributions differ by less than a 
specified tolerance. 

This new update procedure for the displacement thickness was 
originally an ad hoc assumption, but further examination shows 
that it is somewhat analogous to the relationship between changes 
in u Q and 6 given by the von Karman momentum integral 
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( 5 ) 


de 26 + 6* du e 

dx u e dx 



* 

If the shape factor H 5 — is introduced and its streainwise 
variation is neglected, then equation (5) can be approximately 
written as 


* .* T 
d6 » - (H + 2) ~ du + Hdx — ~ 

e PU * 


( 6 ) 


The update procedure on the displacement thickness, including a 
relaxation parameter to, can be written as 


ki+1 


.* i , 

6 + wo 


C e i / 


( 7 ) 


where if to = 1 equation (4) is obtained. With the definition 


Au = u — u 
e u ej e v 


( 8 ) 


and correspondingly 


* *i+l ** 

A5 = 6 - 5 


( 9 ) 


Equation (7) can be rearranged to give the approximate relation 


A6 » - oj — Au 
u e e 


( 10 ) 
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Comparison of equations (6) and (10) shows that they are quite 

similar in rela 4 an increment in u to that in 6 , 

e 

particularly where the skin friction is small. It is recognised 

that this comparison is not exact since the increments in 

* e 
and 6 in the momentum integral equation are changes in the 

streamwise direction whereas those in equation (10) refer to 

changes between successive interaction cycles. Nonetheless, 

there seems to be some merit to this comparison since; for 

example, an alternate 6* update procedure given by 



which is not analogous to the momentum integral equation was 
found to be divergent. Following the same steps which led to 
equation (10), equation (11) can be approximately written as 


AS* a# 


u 


u 


Au 


e 


( 12 ) 


which differs in sign from equation (10). 

Most viscous- inviscid interaction schemes require under- 
relaxation (u> < 1) for convergence, but in the present calcula- 
tions it was found that overrelaxatioa (« > 1) could be used to 
accelerate the convergence of the interaction process. The 
similarity between equations (6) and (10) offers an explanation 
since, even with overrelaxation, it was observed in the present 
calculations that 


(13) 
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w < H + 2 

* 



since H typically varies from 1 to 2.5, while u had to be 
restricted to values less than 2 to avoid divergence. 

RESULTS AND DISCUSSION 

Interaction results obtained with this procedure were 
obtained for the turbulent incompressible flow over ellipsoids 
with 10/1 and 10/3 axis ratios and stings which are 20 percent 
and 33 percent, respectively, of the maximum body radius. A 
small fillet was inserted between the body and sting to avoid 
a slope discontinuity. The initial 5 distribution that was 
used to initialize the interaction procedure was obtained from 
a direct boundary- layer computation using a modified inviscid 
tangential surface velocity as shown in figure 4 for the 10/3 
ellipsoid. This bcundary-layer computation was made with the 
same program that was used in the inverse calculations since, as 
discussed in reference 3, either the direct or inverse option 
can be exercised in this boundary-layer formulation with only a 
simple boundary condition change required. It was necessary to 
flatten the inviscid velocity distribution 3uch that the boundary 
layer did not separate and give a singularity in the initial 
direct boundary- layer computation. Figure 4 also shows, for 
comparison, the final interacted tangential velocity distribution. 

The present calculations are intended only to demonstrate 
this new interaction procedure and thus, to simplify the calcula- 
tions, the interaction was assumed to be negligible upstream of 
the maximum body radius. At this point, the turbulent boundary- 
layer profile was assumed to be given by Coles’ (ref. 13) wall- 
wake formulation with the edge velocity given by the uninteracted 
inviscid computation. Values for the boundary-layer thickness 
and skin-friction coefficient, which are also required in the 
wall-wake formulation, were assumed and are typical of a turbu- 
lent boundary layer at a local Reynolds number of 7.5 x io^ in 
a zero pressure gradient. 
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In the present calculations, 93 grid points were used across 
the boundary layer and 80 ic the streamwise direction. The in- 
viscid calculation was made with 81 and 41 grid points in the 
streamwise and lateral directions, respectively. The downstream 
boundary in the inviscid computation was placed at approximately 
the same location as that used in the boundary-layer computation. 
The velocity potential along this boundary was obtained by quad- 
ratic extrapolation of the upstream values. No problems were 
encountered using this procedure since the boundary was placed 
sufficiently far downstream where the flow velocity had nearly 
recovered to the free-stream value. 

The interacted results for the pressure, displacement body, 

and skin friction are shown in figures 5 and 6 for the 10/1 and 

10/3 ellipsoids, respectively. The Reynolds number used in these 

6 

calculations is 1.53 x IQ based on the body length and free- 
stream velocity. Attached flow was computed for the 10/1 ellip- 
soid as shown in figure 5(b); separated flow, for the 10/3 
ellipsoid, as shown in figure 6(b). The points of separation 
and reattachment are denoted as S and R, respectively, in 
figure 6. In both cases near the upstream boundary the inter- 
acted pressure differs only slightly from that obtained in the 
inviscid calculation and thus it is a reasonable assumption to 
neglect the boundary layer upstream of the maximum body radius. 
Figures 5(a) and 6(a) show typical interaction results, where 
it is seen that the inclusion of the boundary layer results in 
pressure distributions with smaller gradients than the corres- 
ponding inviscid distributions. 

In these calculations a boundary-layer calculation was made 
after every 20 relaxation cycles in the inviscid code RAXB0B. 
Further study needs to be made to determine the optimum number 
of inviscid cycles between boundary-layer computations so as to 
reduce computer time. The 10/1 ellipsoid took about 7 inter- 
actions for convergence and required about 1.5 minutes of CDC 
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CYBER-175 computer time. The 10/3 ellipsoid took about 10 
iterations and 2 minutes of computer time for convergence. 

Conventional viscous-inviscid matching techniques generally 
require, even for mild interaction cases, numerical smoothing 
and underrelaxation after each interaction iteration to prevent 
divergence. In contrast, the present calculations, in which the 
interaction is quite significant, required no smoothing or under- 
relaxation. Smoothing is probably not needed in the present 
interaction procedure since it naturally occurs in the inverse 
boundary-layer calculation when the tangential velocity, u g , is 
computed by integrating the pressure gradient deduced in the cal- 
culation. Numerical integration tends to smooth the results in 
contrast to numerical differentiation, which occurs in direct 
boundary-layer techniques where the imposed pressure gradient is 
computed by numerically differentiating the inviscid pressure 
distribution. 

The convergence of the interaction calculation was acceler- 
ated somewhat as shown in figure 7 by using overrelaxation, 
w = 1.5. In figure 7(a), six streamwise locations along the 
displacement thickness distribution are indicated. In 
figures 7(b) and 7(c) the convergence history of the displace- 
ment thickness at these six locations is shown for both ui = 1 
and w = 1.5. It is seen that the use of overrelaxation has the 
greatest effect in reducing the number of iteration cycles at 
locations c and d, where the interaction is the strongest. Over- 
relaxation was used for only five interaction cycles as it was 
found that an oscillatory pattern in the convergence history 
persisted if overrelaxation was used throughout. Although the 
effect shewn here of overrelaxation is relatively small, it does 
appear that the present interaction procedure is very promising. 
Further study is needed to deduce an optimum value of the over- 
relaxation parameter. 
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CONCLUDING REMARKS 


Converged solutions obtained with this new interaction 
procedure demonstrate that this technique may be quite useful 
for flows with strong interaction including small separated 
regions. The present approach may even be preferable to existing 
interaction schemes for attached flow since the present proce- 
dure, in the calculations reported, required no smoothing and 
over-relaxation was used. This procedure should now be extended 
to compressible flow and comparisons made with experimental data. 
Since the present boundary- layer approach is a finite-difference 
technique, the present interaction procedure should be quite 
useful for testing various turbulence models to improve the 
current capability in predicting the detailed characteristics 
of separated turbulent flow. 
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(a) Pressure distribution and displacement body. 

Figure 5.- Interaction results for 10/1 ellipsoid w/20% sting. 





Figure 5.- Concluded 




(a) Pressure distribution and displacement body. 

Figure 3.- Interaction results for 10/3 ellipsoid with 33% sting. 
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